      IMPLICIT REAL*8 (A-H, L-Z)
      real m1,m2,tb0,mx,mx2,tbx,eccx,t1,t2,Num1,Num2,Num3
      real N1(20000),N2(20000),N3(20000)
      integer I,J,K
      OPEN(14,file='ns.dat',status='unknown')
      open(100,file='ns1.8.dat',status='unknown')
      N1(0) = 0
      do 100 I = 1, 100000
      read(14,*)m1,m2,tb0,mx,mx2,tbx,eccx,t1
      N1(I) = 7.96e-6*m2*m1**-2.7
      tbx = tbx*(1-eccx**2)**1.5
      tb = 0.11633+0.2928*mx-0.0661*mx**2+0.00745*mx**3-0.0003038*mx**4
      if(tbx.ge.tb)then
c      N1(I) = N1(I)+N1(I-1)
c      N1(I)=2*3.43e-7*m2*m1**-2.7


c       f1 = -170.21+203.8*mx-60*mx**2
c       f2 = -1.3+1.6*mx
c       f3 = 1601-1057.5*mx+175*mx**2
c       f4 = 256530-367228.5714*mx+131428.5714*mx**2
c       f5 = 178-120*mx
c       f6 = 7.8046-5.6076*mx+1.9039*mx**2
c      if((mx.le.1.5.and.tbx.ge.0.1.or.mx.gt.1.5
c     & .and.mx.le.1.65.and.tbx.ge.f1.or.mx.gt.1.65
c     &  .and.mx.le.2.5.and.tbx.ge.2.7.or.mx.gt.2.5.and.
c     &  mx.le.3.and.tbx.ge.f2.or.mx.gt.3.and.mx.le.3.2
c     & .and.tbx.ge.f3).and.(mx.le.1.33.and.tbx.le.600
c     &  .or.mx.gt.1.33.and.mx.le.1.4.and.tbx.le.f4
c     & .or.mx.gt.1.4.and.mx.le.1.45.and.tbx.le.f5
c     & .or.mx.gt.1.45.and.mx.le.3.2.and.tbx.le.f6))then
   
   
       f1 = -125.947+94.643*mx-17.516*mx**2
       f2 = -2.9+1.2*mx
       f3 = -11+3*mx
       f8 = -60.06+12.79*mx
       f4 = -3308+620*mx
       f5 = 9007.54-3984.62*mx
       f6 = 13487.54-11153.85*mx+2307.69*mx**2
       f9 = 40.25-12.5*mx
       f7 = 10.312-6.133*mx+2.2064*mx**2
       
       
      if((mx.le.2.42.and.tbx.ge.0.1.or.mx.gt.2.42
     & .and.mx.le.2.7.and.tbx.ge.f1.or.mx.gt.2.7
     &  .and.mx.le.4.and.tbx.ge.1.9.or.mx.gt.4.and.
     &  mx.le.4.5.and.tbx.ge.f2.or.mx.gt.4.5.and.mx.le.5
     & .and.tbx.ge.f3.or.mx.gt.5.and.mx.le.5.35.and.tbx.ge.f8
     & .or.mx.gt.5.35.and.mx.le.5.4.and.
     &  tbx.gt.f4).and.(mx.le.2.11.and.tbx.le.600
     &  .or.mx.gt.2.11.and.mx.le.2.24.and.tbx.le.f5
     & .or.mx.gt.2.24.and.mx.le.2.42.and.tbx.le.f6
     & .or.mx.gt.2.42.and.mx.le.2.5.and.tbx.le.f9.or.
     &  mx.gt.2.5.and.mx.le.5.4.and.tbx.le.f7))then   
c       if(mx.le.6.and.(mx2.ne.1.30.or.m1.gt.9.0))then
       write(100,12) m1,m2,tb0,mx,mx2,tbx,eccx,t1,N1(I)
        endif
100      endif
12    format(8f8.2,1E10.2)
       end
